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Abstract: Bauer, Ligeti and Luke have recently proposed a new method for 
measuring \V u b\ in inclusive semileptonic B decays, using a cut ^fq 1 > Mb — Mb on 
the lepton invariant mass to discriminate against b — > c transitions. We investigate 
the structure of the heavy-quark expansion for this case and show that to all 
orders the magnitude of the leading perturbative and nonperturbative corrections 
is controlled by a hadronic scale /i c ^ m c depending on the minimal value of q 2 . 
These corrections can be analyzed using a modified version of the heavy-quark 
expansion ("hybrid expansion"). We find that the theoretical uncertainty in the 
extraction of \V u b\ is a factor 2.5 larger than previously estimated, which allows for 
a determination with 10% accuracy. 
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1. Introduction 

A precise knowledge of the magnitude of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element V u b is of great importance to the study of CP violation at the B fac- 
tories. Whereas the main target of CP-asymmetry measurements is to fix the angles 
of the unitarity triangle, the length of the sides of this triangle are determined by 
the smallest CKM elements, \V u b\ and \Vtd\- The measurement of \V u b\, which neither 
involves CP violation nor rare loop processes, appears to be the simplest step in the 
overall determination of the unitarity triangle. Yet, at present this measurement is 
limited by uncomfortably large theoretical uncertainties. 

\V u b\ can be determined most directly from semileptonic B decays into charmless 
final states. The theoretical interpretation of exclusive decays such as B — > n I v or 
B — > pi v is limited by the necessity to predict the B — > n or B — > p transition form 
factors, which parameterize the complicated hadronic interactions relevant to these 
decays. The inclusive decays B — > Xlv admit a cleaner theoretical analysis based 
on a heavy-quark expansion [1, 2, 3]. However, an obstacle is that experimentally 
it is necessary to impose restrictive cuts to suppress the background from B — > 
X c lu decays (i.e., decays into final states with charm). Accounting for such cuts 
theoretically is difficult and usually introduces significant uncertainties. 

The first determination of \V u b\ from inclusive decays was based on a measure- 
ment of the charged-lepton energy spectrum close to the endpoint region, which is 
kinematically forbidden for decays with a charm hadron in the final state. This re- 
striction eliminates about 90% of the B — > X u I v signal, and it is difficult to calculate 
reliably the small fraction of the remaining events. The reason is that in this por- 
tion of phase space the conventional heavy-quark expansion breaks down and must 
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be replaced by a twist expansion, in which an infinite tower of local operators is 
resummed into a shape function describing the light-cone momentum distribution of 
the b quark inside a B meson [4, 5, 6]. A promising strategy is to infer the fraction 
of semileptonic decays in the endpoint region from a study of the photon spectrum 
in B — > X s 7 decays, since the leading nonperturbative effects in the two decays are 
described by the same shape function [6, 7, 8, 9]. Recently, it has been emphasized 
that a cut Mx < Md on the hadronic invariant mass would provide for a better 
discrimination between B — > X u lv and B — > X c lv decays [10, 11, 12]. Such a 
cut eliminates the charm background, while affecting only about 20% of the signal 
events. Despite of this advantage, however, it is difficult to calculate precisely what 
this fraction is [13]. 

Bauer, Ligeti and Luke (BLL) have pointed out that the situation may be better 
if a discrimination based on a cut on the lepton invariant mass ^Jq 2 is employed [14]. 
Requiring q 2 > (Mb — Mb) 2 eliminates the charm background, while containing 
about 20% of the signal events. This fraction is much less than in the case of a cut 
on hadronic invariant mass, but the inclusive rate with a lepton invariant mass cut 
offers the advantage of being calculable using a conventional heavy-quark expansion, 
without a need to resum an infinite series of nonperturbative corrections. BLL find 
that the fraction of B — > X u lv events with q 2 > q 2 , where ql > (M B — M D ) 2 is 
required to eliminate the charm background, is given by 

F^l) = (i + % 2 )(i - & + x(?o) - (% 2 - 2g 6 + g 8 ) ^ + • • • , (Li) 

where q 2 . = q 2 jm\, and the dots represent corrections of higher order in the expansion 
in powers of a s (m&) and A/m b (with A a characteristic hadronic scale). The function 
X(q 2 ) can be obtained from results presented in [15]. The nonperturbative parameter 
A 2 = \{M B * —M B ) ~ 0.12 GeV 2 was introduced in [16]. As it stands, eq. (1.1) seems 
to provide a solid theoretical basis for a systematic analysis of the partial decay rate 
in an expansion in logarithms and powers of Kjm^. (This leaves aside the fact that 
the lepton invariant mass cut eliminates about 80% of all B — > X u I v events, and 
therefore one may worry that violations of quark-hadron duality may be larger than 
for the total inclusive semileptonic rate.) 

The purpose of this work is to analyze the structure of the result (1.1) in more 
detail. We show that the relevant mass scale /i c controlling the size of corrections 
in the heavy-quark expansion is less than or of order the charm-quark mass, rather 
than the heavier 6-quark mass. This is so because the largest values of the hadronic 
invariant mass and energy accessible are of order the charm mass or less. Although 
the ratio m c /mb is usually taken as a constant when discussing the heavy-quark limit, 
it is well known that the convergence of heavy-quark expansions at the charm scale 
can be poor. This makes our observation relevant. We suggest that an appropriate 
framework in which to investigate the leading corrections to the heavy-quark limit 
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in the present case is a modified version of the heavy-quark expansion to which 
we refer as "hybrid expansion". The idea is that in the kinematic region where 
M B 3> E x ~ M x 3> A one can perform a two-step expansion in the ratios E X /M B 
and A/ Ex- A similar strategy has been used in applications of the heavy-quark 
effective theory (HQET) to resum the so-called "hybrid logarithms" a s ln(m c /mb) 
[17] arising in current-induced b — > c transitions using renormalization-group (RG) 
equations [18, 19]. In the context of inclusive decays, an approach similar in spirit 
to our proposal was suggested first by Mannel [20]. We use the hybrid expansion 
to obtain a RG-improved expression for the perturbative corrections to the quantity 
i*X?o) at next-to- leading order (NLO), as well as to estimate the size of higher-order 
power corrections omitted in (1.1). 

2. Structure of the hadronic tensor 

The strong-interaction dynamics relevant to inclusive B decays is encoded in a 
hadronic tensor defined as the forward matrix element of the time-ordered prod- 
uct of two weak currents between S-meson states. Although the variable of prime 
interest to our discussion is the lepton invariant mass, the physics of the hadronic 
tensor is most naturally described in terms of the hadronic invariant mass and en- 
ergy, M x and E x . These variables are related by q 2 = M| — 2M B E X + M x , and 
thus the restriction q 2 > q^ implies 

Mx<^x< Ml+ 2 ^~ 9 ° 2 , Ml<M 2 x <(M B -^) 2 . (2.1) 

With the optimal choice q 2 = (M B - M D ) 2 this gives M, < M x < M D and M x < 
E x < M D - \{Mj) - M X )/M B . Both variables vary between M n and M D . If the 
cutoff q 2 is chosen higher, as may be required for experimental reasons, their maximal 
values become less than M D . 

The corresponding variables entering a partonic description of inclusive decay 
rates are the parton invariant mass and energy, ^fp 1 and v • p, where p is related to 
the lepton momentum q by p = m^v — q, and v is the velocity of the B meson. The 
phase space for the dimensionless variables p 2 = p 2 /m 2 and z — 2v • p/rrib is 

2V^< z < 1 -q 2 +p 2 , <p 2 < (1 - ygf) 2 . (2.2) 

For the purpose of our discussion here is the pole mass of the b quark. Alternative 
mass definitions will be discussed in more detail later. The optimal value of q 2 is 
q 2 = (M B — M B ) 2 /m 2 ~ (1 — — ) 2 . In this case the largest values of the parton 
variables are z max = 2 ^ and p 2 nax = (^) 2 - Without a restriction on q 2 , the phase 
space for these variables would be such that z max = 2 and p„ iax = 1. In other words, 
the lepton invariant mass cut restricts both variables to a region where they are 
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parametrically suppressed, such that v ■ p and \fp 1 are at most of order m c . Besides 
nib, the charm-quark mass is then a relevant scale that determines the magnitude of 
strong-interaction effects in the heavy-quark expansion. 

To make the parametric suppression noted above more explicit, we introduce a 
characteristic scale /i c and an associated small expansion parameter e by 



He 



2m h 



0(m r 



(2.3) 



and rewrite the fraction F(q^) of B — > X u I v events with q 2 > q 2 as 



F(e) = 16e 3 (l - e) 



C (6)-|c 2 (e)^ + O[(A//, c ) 3 ] 

z H'c 



(2.4) 



where C n (e) = 1 + 0(e, a s ) are short-distance coefficients. From (1.1) we find 

-)+e(-81n2 e + ^)+0(e 2 ) + 0(a 2 ) , 

(2.5) 



-61n2e 

2 3 



^/a l-8e 2 + 8e 3 . 
C 2 (e) = + 0(a s ) . 



1-e 

The definition of e and /i c , as well as the explicit form of the coefficients C n (e), 
depend on the definition of the heavy-quark mass. The above result for Co(e) refers 
to the pole mass. Later we will introduce a more suitable mass definition and give 
an exact NLO expression for Co(e) including the higher-order terms in e. 

It is apparent from (2.4) that the leading 
power correction proportional to A 2 is of order 
(A//x c ) 2 . It is not difficult to see that, in the 
presence of a lepton invariant mass cut, also in 
higher orders the power corrections scale like 
(A//x c ) n . For simplicity of the argument we 
work to leading order in a s , where only the tree 
diagram shown in Figure 1 contributes to the 
hadronic tensor. In momentum space, the u- 
quark propagator gives a contribution 



V 




p+k 




Figure 1: Tree- level contribution 
to the hadronic tensor. Dashed lines 
represent the weak currents. 



(p + ky 



p 2 + 2p ■ k + k 2 



(2.6) 



where k = 0(A) is the residual momentum of 
the heavy quark inside the B meson [21]. Roughly speaking, the heavy-quark ex- 
pansion is obtained by replacing the residual momentum with a covariant derivative, 
k^ 1 — ■> iD^, thereby introducing the soft interactions of the w-quark jet with the back- 
ground field of the light degrees of freedom in the B meson. Let us discuss how the 
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three terms in the denominator of the propagator scale in the different kinematic 
regions relevant to the determination of \V ub \. Whereas the k 2 ~ A 2 term is always 
suppressed, the relative magnitude of the two other terms depends on the kinematic 
region considered. For the regions of the large charged-lepton energy or low hadronic 
invariant mass, the terms p 2 ~ maxfAroj, m 2 } and p- k ~ Amj are of the same mag- 
nitude, and it is thus necessary to resum terms of the form (p ■ k/p 2 ) n to all orders in 
the heavy-quark expansion. This leads to a twist expansion, where these terms are 
absorbed into a nonperturbative shape function [4, 5, 6]. In contrast, for large lepton 
invariant mass, p ■ k ~ A m c is parametrically suppressed with respect to p 2 ~ m 2 by 
a power of A/m c . 

These general observations can also be derived from the explicit expression for 
the differential decay rate expressed in terms of the variables z and p 2 , normalized 
to the total decay rate. At NLO in the heavy-quark expansion we find 

1 iM 2z 2 (3-2z)5(p 2 ) + ^E(z,p 2 ) 



T dzdp 2 



3tt 



- KP 2 ) 



- np 2 ) 



Ai 
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mi 
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mi 
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mi 



mi 



Ai 



- *"(P 2 ) — (3 — 2z) — ^ + 
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(2.7) 



where Ai is a nonperturbative parameter related to the average kinetic energy of the 
b quark inside the B meson [16], and the function E(z,p 2 ) gives the perturbative 
correction calculated in [13], which also includes the correction to the total decay 
rate appearing in the denominator on the left-hand side. The power corrections in 
(2.7) have been derived using the results of [3, 22]. In the kinematic region where 
z = 0(e) and p 2 = 0(e 2 ), it is instructive to change variables from (z,p 2 ) to (z, £)> 
where £ = Ap 2 /z 2 = p 2 / (v ■ p) 2 G [0, 1] is related to the parton velocity in the B rest 
frame. The kinematic range for these variables is 



o < e < i, 



< z < - (1 - yjl - 2e0 



2e + e 2 e + 0(e 3 ) 



(2.8) 



The variable £ is of order unity irrespective of the lepton invariant mass cut. Thus, 
after the transformation the only parametrically small quantity is z — 0(e). In terms 
of the new variables, the double-differential decay rate turns into an expansion in 
powers of A/(zm&), and the perturbative corrections contain single logarithms of z. 
Integrating the double-differential rate over z, and keeping only terms of leading 
order in e, we obtain 
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(2.9) 



The integral of this expression over £ reproduces the leading terms in e in (2.4). The 
result for the power corrections shows that indeed /i c = 0(m c ) is the characteristic 
scale of the hybrid expansion. Perturbative logarithms of fi c / m b appear because 
the intrinsic scale of the hadronic tensor is smaller than the mass of the external b 
quarks, suggesting that the appropriate scale to evaluate the running coupling a s 
is significantly less than m b . This is in accordance with the observation [23] that 
the physical scale derived using the Brodsky-Lepage-Mackenzie (BLM) scale-setting 
prescription [24] strongly decreases with increasing q 2 . We will come back to the 
question of scale setting in the next section. 

For later purposes, it will be useful to have 
yet another way of reproducing the result (2.4). 
To this end, we represent the fraction F(q^) as a 
contour integral in the complex p 2 plane. Such 
a representation exists because the hadronic ten- 
sor is an analytic function in the complex plane 
apart from discontinuities located on the posi- 
tive real p 2 axis, as illustrated in Figure 2. We 
write 




n<$) = ^ f dp 2 T(p 2 ,q 2 ), (2.10) 



Figure 2: Contour representation 
of F(<7o) in the complex p 2 plane. 



|P 2 |=Pmax 



where p^ ax = (1 — y<?o) 2 , and the correlator 
T{p 2 ,qo) can be obtained using dispersion re- 
lations [13]. Eliminating q 2 in favor of e, and parameterizing p 2 = p^ ax e lv on the 
contour of integration, the result can be written in the form 



2n 



F( e ) = 16e 3 (l - e) J ^(e'V) , 




dip 



(2.11) 



where the function t(e t(fi , e) contains the corrections to the heavy-quark limit. We 
obtain 



f (e* e) = C (e) + |i Y(e* e) - § C 2 (e) *§ + 



H'c r c 



-vp 



(2.12) 
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The function Y(e i(f , e) with 

Vfe e) = T * - — + ( T - - *) „ln — - 2 In' — + 0(e) (2.13) 

and u = \/l — x is defined such that its contour integral vanishes. The coefficients 
Di(e) and £i(e) are equal to 1 in the limit e — > 0. Their explicit expressions are 

(i - e)( i - vr^)' (i-r + r 2 ) + 0(Q ->- 

__A__( 1+ 3 £ _ 4f2 ) +OK)I 

d-Od-Vi^' ^-T^ ^'- (2 ' 14) 

Equations (2.9) and (2.12) allow us to study the behavior of the leading cor- 
rections in the heavy-quark expansion in more detail. In particular, we will utilize 
them to estimate unknown, higher-order power corrections. We will, however, first 
investigate the perturbative corrections in more detail, using the hybrid expansion 
as a tool to perform a systematic RG improvement of the one-loop expressions in 
(1.1) and (2.5). 



Di(e) = 
D 2 (e) = 



3. RG improvement and definition of the heavy-quark mass 

Based on the observation that the event fraction F{q^) in (1.1) receives very small 
corrections of order (3qOl 2 s , where f3 is the first coefficient of the /3-function, BLL have 
argued that the perturbative uncertainty in their prediction is negligible [14]. The 
purpose of this section is to critically reanalyze the perturbative uncertainty in the 
calculation of this quantity. We first note that in the present case it is misleading to 
associate the size of (3 a 2 corrections with a physical BLM scale in the process. The 
total semileptonic rate and the lepton invariant mass spectrum receive very large 
corrections of order (3oa 2 s) corresponding to very low BLM scales [23]. The fraction 
F(e) is defined as the ratio of the partially integrated lepton spectrum and the total 
rate. If both of these quantities have low physical scales, the same must be true for 
their ratio. In our opinion the small fact 2 , correction observed in [14] is thus due to 
an accidental cancellation and does not bare any physical significance. 

Here we follow a different strategy to estimate the potential importance of higher- 
order effects. We have argued that a useful framework in which to analyze the fraction 
F(e) is provided by a hybrid expansion, in which the physics associated with the 
three mass scales mj ^> /i c ^> A is disentangled. Since the intermediate scale fi c is 
only about 1 GeV or less, and since the running of the strong coupling in the region 
between m& and 1 GeV is significant, we expect important higher-order perturbative 
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corrections resolving the scale ambiguity. At one-loop order, this is indicated by the 
presence of logarithms of e in (2.5). At any given order in an expansion in powers 
of e the contributions associated with the two couplings « s (m(,) and a s (fj, c ) can be 
separated by solving RG equations in the hybrid expansion at NLO. The residual 
scale ambiguity left after RG improvement provides an estimate of the perturbative 
uncertainty in the result. Unfortunately, to perform this program one must compute, 
at every order in e, the two-loop anomalous dimensions of a new tower of higher- 
dimensional operators. At present, these anomalous dimensions are known only for 
the operators entering at the leading order in e, although partial results exist for the 
operators relevant to the 0(e) terms [25]. 

We now discuss the RG improvement at leading order in e in detail. The first 
step in the construction of the hybrid expansion is to expand the weak currents in 
the definition of the hadronic tensor in terms of operators of the HQET, with the 
result [21] 

^(1 - 75 ) d (^) q^(l - 75) h v + C 2 (—) qv^l + 75) K + 0(l/m b ) , (3.1) 

where is the _B-meson velocity, h v are the velocity-dependent fields of the HQET, 
and [i is the scale at which the operators are renormalized. The RG-improved ex- 
pressions for the Wilson coefficients Ci(rrib/fi) are known at NLO. The terms of order 
l/m b in (3.1) would contribute at order e in the hybrid expansion and can be ne- 
glected for the discussion of the leading terms. It is important that all dependence 
on the 6-quark mass is explicit in (3.1). We now insert this result into the time- 
ordered product of currents in the hadronic tensor and perform an operator product 
expansion of the current product. This is an expansion in logarithms and powers of 
A/ fj, c . The scale nib does not appear in the matrix elements of the hybrid expansion. 
We find that at leading order in e only the HQET current product proportional to 
C\{rribl \i) contributes. Using the known NLO expression for the coefficient C\{mbj '/i) 
[26, 27], we obtain 



C (e) 



Here (3 



a 



a 



25 



a s (m b ) ) 



a s (m b ) ) 
and Z hl = 



1 + 



2a s (m b ) 
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+ 0(e) 



+ 0(e). (3.2) 



— !|jr are perturbative coefficients evaluated for rif — 4 



3 m 7500 

light quark flavors, as is appropriate for a scale of order m c . Note that the NLO 
corrections proportional to the coupling a s (mb) contain the corrections to the total 
semileptonic rate. The above result for Co(e) gives the RG-improved form of the 
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leading term in e in the one-loop expression in (2.5). This result is formally scale 
independent at NLO. The renormalization scale [i should be chosen of order \i c in 
order to avoid large logarithms in the hybrid expansion. (The inappropriate choice 
fi = nib would reproduce the one-loop result with a s evaluated at nif,.) 

The dashed line in Figure 3 shows 
the perturbative prediction for the frac- 
tion F(e) at leading order in e, as a func- 
tion of the renormalization scale. Here 
and below we use the two-loop expression 
for the running coupling constant nor- 
malized such that a s (Mz) = 0.118. In 
accordance with relation (3.4) below we 
use nib = 5.0 GeV for the pole mass, not- 
ing that this value (and thus the normal- 
ization of the dashed curve) has a large 
uncertainty. Two observations are im- 
portant. First, the scale dependence of 
the dashed curve is significant, and hence 
the result obtained with an appropriate 
choice of scale /i m fi c is much lower than 
that obtained with the naive choice fi = nib- Secondly, perturbation theory in the 
on-shell scheme breaks down at a scale not much less than the appropriate scale 
jj, m /j, c m 1 GeV. We conclude that in the on-shell scheme there is a large pertur- 
bative uncertainty in the calculation of the coefficient Co(e), which is not apparent 
from the naive one- loop result. This conclusion is in contrast with the assumption 
made by BLL, that the perturbative uncertainty is negligible [14]. 

The breakdown of perturbation theory at a scale of order /i c can be traced back 
to the large coefficient of the NLO correction proportional to a s (fi) in (3.2). We will 
now show that the size of this coefficient can be reduced significantly by adopting 
a more appropriate definition of the heavy-quark mass. So far we have worked in 
the on-shell scheme, where the mass is defined as the pole in the renormalized quark 
propagator, rrib = itlI o1b . Since the pole mass is affected by IR renormalon ambiguities 
[28, 29], it is better to eliminate it from the final expressions for inclusive decay rates. 



Figure 3: Scale dependence of the pertur- 
bative result for the fraction F(e) to leading 
order in the hybrid expansion, evaluated for 
q 2 = {Mb — Mo) 2 - The two curves refer to 
different mass renormalization schemes. 



If a new mass definition m' h is introduced via m£ ole 
that 



Z m m' b , it follows from (2.4) 



C' (e' 



a 



pole / 



Z m -\ (1-2(0(3 -460 



e' 



1-e' 



(3.3) 



where the prime on e indicates that this parameter in sensitive to the definition of 
m b . We observe that a multiplicative redefinition of m& with Z m (a s ) = 1 + 0(a s ), 
such as the relation between the pole mass and the running mass defined in the MS 
scheme, is not appropriate in our case, since it would lead to a contribution to C (e') 
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that is enhanced by a factor of a s /e'. This would upset the power counting in the 
hybrid expansion. We suggest instead to work with a short-distance mass subtracted 
at a scale /if of order /i c , which is the natural scale of our problem. 

It is well known that the convergence of the perturbative series for near on-shell 
problems in heavy-quark physics can be largely improved by introducing low-scale 
subtracted quark masses, which have the generic property that they differ from the 
pole mass by an amount proportional to the subtraction scale fif. Several such mass 
definitions exist and have been applied to various processes [30, 31, 32, 33]. To 
illustrate the point, we use the potential-subtracted (PS) mass m PS (//y) introduced 
by Beneke [30] and evaluate it at the scale /i/ = /i^ 8 . 1 At NLO, the relation between 
the pole mass and the PS mass reads 



m 



pole 



PS 



3tt 



1 + 



2tt 




A) In 



-4 



+ 



, (3.4) 



which is formally independent of the scale fi at which the coupling is renormalized. 
Note that the difference between the two mass definitions is a perturbative series 
multiplying the scale yU PS = e PS m PS . At NLO in a s , it then follows from (3.3) that 



C PS (e PS ) = C polc (e PS ) 



1 + 



4a s (/i) /i PS (l-2e PS )(3-4e PS ) 



3tt e p Sm PS 



e/ PS^ 



1 + 4^) +0(ePS) 



7T 



C PS 



(3.5) 



From now on we will use the PS mass m b = m PS (/i PS ) in all our equations and omit 
the label "PS" on the quantities m b , fi c and e. Using (3.2) and adding the extra 
contribution proportional to a s (fj,), we obtain 



C PS i 



e ~ 



a s (m b ) 



12 

25 



0.71 



a s {m b/ 



7T 



+ ^Mfo.63 + 21n^ 



+ 0(e). (3.6) 



The introduction of the PS mass has much reduced the size of the NLO correction. 
The result for the leading contribution to C*o(e) in the PS scheme is shown by the 
solid line in Figure 3. It exhibits a better stability that in the on-shell scheme, and 
it is stable down to lower values of the renormalization scale. 

The value of the PS mass at the scale /i 2 = 2 GeV has been determined from a 
sum-rule analysis of the bb production cross section near threshold, with the result 
m PS (2 GeV) = (4.59 ± 0.08) GeV (corresponding to m b (m b ) = (4.25 ± 0.08) GeV in 
the MS scheme) [34]. At NLO, we can use relation (3.4) to convert this into a value 



1 We could instead evaluate the PS mass at a scale /!/ = £ yu c with £ = 0(1), however this would 
lead to more complicated expressions. Varying £ between 1 and 2 leads to a variation of the results 
by an amount similar to the perturbative uncertainty estimated later in this section. 
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of the PS mass at the scale /i c . This gives the implicit equation 



ml s (ii c ) = m^ s (/i 2 ) +/i 



PS/ 



4ct s (/x 2 
3tt 




25 a s (n 2 ) ln A*2 
(i 11 



7T 



(3.7) 



from which we determine the scale \i c and then the mass m b = m^ s (/i c ). For example, 
we find m h ~ 4.73 GeV, /i c ~ 1.13 GeV, e ~ 0.24 for g 2 = (m b - M D ) 2 , and 
m b ~ 4.79 GeV, /i c ~ 0.83 GeV, e ~ 0.17 for gg = 15 GeV 2 . 

In (3.6) we have obtained a RG-improved expression for the short- distance co- 
efficient at leading order in e. It is at present not possible to extend this analysis to 
higher orders in the hybrid expansion, since the corresponding two-loop anomalous 
dimensions of higher-dimensional operators are unknown. However, since in the PS 
scheme the leading term in e gives the dominant contribution to Co(e), we expect 
that the unresolved scale ambiguity in the higher-order terms does not introduce a 
large uncertainty. Our final expression for the short-distance coefficient at NLO is 



C n ps 



a s (m b ) J 



Q 71 a s (m b ) a s (n) 



71 



TV 



0.63 + 21n-^) +—G(e) 



(3.8) 

where the scale in a s in the 0(e) term remains undetermined. The exact result for 
the function G(e) in the PS scheme is 



G(e) 



1 



12(1 



1 



4 



= — In 2e 



6(1 -e)e 3 
L 2 (2e) 
95 



1 89 9 , 

12-12 e + 21e +4 ° e 



64e 4 + 



1 + e + -e 2 + 2e l 



ln2e 



L 2 (l 



2f » + T 



18 



32 



1337\ 



+ H-i5"' 2f+ wJ +0(f) 



(3.9) 



The left-hand plot in Figure 4 shows Co(e) for the optimal choice g 2 = (Mb — Mp) 2 as 
a function of the renormalization scale. The width of the band reflects the sensitivity 
of the result to the value of the coupling d s associated with the 0(e) terms in (3.8), 
which we vary between a s (m h ) and a s (/i). To estimate the residual scale dependence 
we vary \i between the values \i c and 2/i c . For lower values the perturbative expansion 
diverges, since the running coupling a s (ji) strongly increases below // 1 GeV. 
For comparison, we mention that the naive perturbative analysis with fixed scale 
fi = rrib adopted in [14] would give the much smaller value Co(e) ~ 1.16 at minimal 
gg — 11.6 GeV 2 . The fact that we find larger QCD corrections will have important 
implications for the extraction of \V u b\. The right-hand plot in Figure 4 shows C (e) 
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yu [GeV] q 2 [GeV 2 ] 

Figure 4: Left: Scale dependence of the coefficient Cq (e) for q 2 . = (Mb — Mr,) 2 . 
The width of the band indicates the scale ambiguity of higher-order terms in the hybrid 
expansion. The vertical lines show the window between fj, = fj, c and fi = 2[i c used to 
estimate the scale dependence. Right: Coefficient Co(e) including theoretical uncertainties, 
as a function of the lepton invariant mass cut (see text). 

as a function of q 2 . The width of the band represents the total scale uncertainty, 
estimated by variation of \i and a s as described above. 

An independent way to estimate the uncertainty in the value of the perturba- 
tive coefficient Co(e) is based on the contour representation (2.11) for the fraction 
F(e). As we have just discussed, the value of C (e) depends on the definition of the 
heavy-quark mass. However, the variation of the 0(a s ) correction in (2.12) along 
the circle in the complex momentum plane is independent of mass redefinitions. We 
may thus take the ^-variation of the one-loop correction, given by a s (/i c )/3iT times 
the variation of the real part of the function Y(e lip ,e) in (2.13), as a typical size of 
an 0(a s ) correction in the problem at hand. For an asymptotic series, the value 
of that correction provides an estimate for the magnitude of unknown higher-order 
corrections. The dashed lines in the right-hand plot in Figure 4 show this variation 
as an error band applied to the central values of C (e). This independent evalua- 
tion of higher-order effects is in good agreement with our previous estimate of the 
perturbative uncertainty, giving us confidence that this estimate is a realistic one. 

4. Higher-order power corrections 

Uncertainties enter the theoretical prediction for the fraction F(e) also at the level 
of power correction. First, there are unknown 0(a s ) corrections to the Wilson coef- 
ficient 02(e) multiplying the term proportional to A2/ 1 \j? c in (2.4). To estimate their 
effect, we replace the bracket [. . .] in this equation with C (e) [1 — § C^e) X2/ fJ^ + ■ • 
which amounts to multiplying the tree- level coefficient C^e) in the original expres- 
sion with C (e). At q = (M B — M D ) 2 , the difference is a 5% effect. Potentially 
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more important are higher-order power corrections scaling as (A//i c ) 3 . The opera- 
tor matrix elements contributing at third order in the heavy-quark expansion can be 
identified [35, 36], but little is known about their actual size. Naive dimensional anal- 
ysis suggests that, with a typical hadronic scale A « 0.5 GeV, a third-order power 
correction could be of order (A//i c ) 3 ~ 0.09 for g 2 = (M B -M D ) 2 and (A//i c ) 3 ~ 0.22 
for g 2 = 15 GeV 2 , but clearly these are rough estimates which must be taken with 
caution. 

We will attempt to extract as much 
information about power corrections as 
possible from the formulae derived in Sec- 
tion 2 for the B — > XI v decay rate and 
spectra in the presence of a lepton invari- 
ant mass cut. We start with the quan- 
tity F(e) itself, which as shown in (2.4) 
receives a moderate second-order power 
correction proportional to the parameter 
A 2 ~ 0.12 GeV 2 . The dashed line in Fig- 
ure 5 shows an estimate of the unknown 
(A/ He) 3 correction, obtained by raising 
this second-order term to the power 3/2. 
To address the question to what extent 
this is a conservative estimate of a "generic" higher-order correction, we focus on 
the differential spectrum in (2.9) and on the contour representation in (2.11). The 
function F(e) is obtained from these results by performing integrals over £ or over 
the contour in the complex plane, respectively. However, the differential distribu- 
tions contain additional information about power corrections, which is not seen after 
the integrations are performed. We first discuss the case of the contour integral in 
(2.11), taking the point of view that except for the region of small (p the magnitude of 
t(e lip , e) can be used to estimate the "generic" size of corrections to the heavy-quark 
limit. This is so because on any point on the circle far away from the real, positive 
p 2 axis, the function t(e lip , e) admits an operator product expansion in a series of 
local operators, whose contributions scale like powers of A//i c . The average over the 
circle determines the corrections to the function F^q^). However, the finer details of 
the distribution on the circle become relevant, e.g., when in a real experiment events 
with different hadronic masses and energies are weighted by different efficiencies. In 
other words, the various terms proportional to A1//1 2 and A2//i 2 in (2.12) are as valid 
as estimate of a second-order power correction as is the \il term in (2.4). Specifi- 
cally, we calculate the average value of the modulus of the power corrections in (2.12) 
on the circle in the complex plane, and then we raise this number to the power 3/2 to 
obtain an estimate of a "generic" (A//i c ) 3 correction. For our numerical analysis we 
use Ai = — (0.30±0.15) GeV 2 , which is in the ball park of recent determinations [37]. 




11.5 12 12.5 13 13.5 14 14.5 15 



qo 2 [GeV 2 ] 

Figure 5: Estimates of third-order power 
corrections to the function F(e), in units of 
the tree- level contribution (see text). 
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The result is shown by the dark band in Figure 5, whose width reflects the sensitivity 
to the value of Ai. If we were to consider larger values of |Ai| the upper limit of the 
band would increase. Another estimate of power corrections can be obtained from 
the coefficients of the various 5-function terms in (2.9). If we take one third of the 
geometric average of the three coefficients, and raise the result to the power 3/2, we 
obtain the light band shown in the figure. 

The above analysis shows that, as expected, there is a large uncertainty in the 
estimate of higher-order power corrections. We do not claim that these corrections 
are likely to be as large as indicated by the upper limit of the light band in Figure 5, 
however, as we have shown this would indeed be possible without introducing any 
unnaturally large coefficients or parameters. Keeping this caveat in mind, we will 
from now on use the upper limit of the dark band in the figure as our estimate of 
third-order power corrections. Numerically, this estimate is close to the one obtained 
by BLL [14]. 



5. Phenomenological implications and summary 

The proposal of BLL is to use the theoretical calculation of the fraction F(q%) to 
obtain a model-independent determination of the CKM matrix element jV^I with 
controlled and small theoretical uncertainty [14]. To this end, one uses the relation 
Br(B — > X u I u)\ q 2 >q 2 = F(qQ)T(B — > X u Iu)tb, where tb is the 5-meson lifetime, 
and T(B — > X u lv) is the total semileptonic decay rate into charmless final states. 
This rate can be calculated with high accuracy in terms of a low-scale subtracted 
6-quark mass, including perturbative corrections of order [a s {m h )] 2 [38] and power 
corrections of order (A/m&) 2 . For our purposes, we use the PS mass defined at the 
scale /i 2 = 2 GeV. Then the expression for the total rate is [33] 

T(B -> X u l V) = IQ2TT^ ^ + " pcrt + Vwer) , (5.1) 

where 5 pcrt ~ 0.04 at two-loop order, and 5 pow cr = Al 2 "^ A2 ^ —0.03. The small 
uncertainties in these two quantities are negligible for our numerical analysis below. 
Using these results, we obtain the master formula 

-Br(B^Xlv)\ q2>q2 l.6ps nl 2 



V ub \ = 2.96 x 10 



-3 



10- 3 F'(g 2 ) t b 
where all theoretical uncertainties are contained in the function 



(5.2) 



^M^W'- (5 - 3) 

The mass dependence due to factor [m^ s (/i2)] 5 from the total decay rate is positively 
correlated with the mass dependence of the function F^q 2 ,). As a result, our predic- 
tions for the function F'(q 2 ) become extremely sensitive to the value of the 6-quark 
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F'{ql) 


5m b a 2 a s {A//2 c ) 2 (A/> c ) 3 


{M B - M D ) 2 
13 GeV 2 
15 GeV 2 


0.204 ±0.040 
0.151 ±0.036 
0.090 ±0.032 


16.7% 6.0% 3.0% 8.2% 
18.5% 7.2% 4.7% 12.5% 
22.2% 10.0% 9.4% 23.8% 



Table 1: Predictions for F'{qq) for three values of the lepton invariant mass cut. The 
right-hand portion contains the relative theoretical uncertainties due to the uncertainty in 
the value of the 6-quark mass (assuming bm^ = 80MeV [34]), higher-order perturbative 
corrections (including scale dependence), perturbative corrections to second-order power 
corrections, and higher-order power corrections. 



mass. For practical purposes, this dependence can be parameterized as 



where 



w 9x On - ( M B ~ M D f , . 



lGeV 

In Table 1, we show our final results for the quantity F'{ql) and its theoretical 
uncertainties (as estimated above) for some representative values of q 2 . For compar- 
ison, we note that BLL obtained the values F'{{M B — M D ) 2 ) = 0.169 ± 0.016 and 
F'(15GeV 2 ) = 0.061 ± 0.013, where the dominant theoretical error was assumed to 
be due to higher-order power corrections. Our central values are significantly higher 
because of the larger perturbative correction obtained after RG improvement. Note 
that our error estimates are about 2.5 times as large as those quoted by BLL. The 
difference between the central values of the two calculations is about la of our errors, 
and about 2a of their errors. 

In Figure 6 we show a graphical representation of the fraction F'{q^) and its total 
theoretical uncertainty. This result, together with the master formula (5.2), provides 
the theoretical basis for the determination of \V u b\. The right-hand plot in the figure 
shows the fractional theoretical uncertainty in the result for \ V u b\. Although our error 
estimates are more pessimistic than those presented by BLL, we still conclude that 
their method provides a very promising route for a precise determination of |V^|. For 
a realistic cut on the lepton invariant mass in the vicinity of gj — 12.5 GeV 2 , which 
is about 1 GeV 2 above the optimal value, the theoretical uncertainty in |T4&| is close 
to 10%. A determination with such an accuracy would be a significant improvement 
with respect to the present knowledge of this important parameter. We believe it 
would also be more reliable than a future determination obtained by combining the 
partial decay rates in the endpoint regions of B — > X s 7 and B — > X u I v decay spectra 
[6, 7, 8, 9], which is limited by uncontrollable power corrections of first order in A/mj, 
that violate the factorization of soft and collinear singularities. 
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11.5 12 12.5 13 13.5 14 14.5 15 

2 



qo z [GeV z ] 



11.5 12 12.5 13 13.5 14 14.5 15 
2 rv^T/2n 



qo A [GeV] 



Figure 6: Left: Theoretical prediction for the quantity F'{c^). Right: Relative theoretical 
uncertainty in the extraction of |V^,|. 



According to Table 1, the dominant sources of theoretical uncertainty in the 
extraction of \V U }\ are associated with the sensitivity to the value of the 6-quark 
mass and with unknown, higher-order power corrections. Whereas it is not obvious 
how one should obtain a reliable value for the power corrections, the precision in the 
value of the 6-quark mass can presumably be improved by reducing the theoretical 
uncertainties in the analysis of (bb) bound states. 2 In addition, it would be possible to 
reduce the perturbative uncertainty in the calculation in two ways, by calculating the 
exact 0(a 2 ) corrections to the fraction F(q^) (the two-loop corrections to the total 
decay rate are known [38]), and by computing the two- loop anomalous dimensions 
of the operators contributing at 0(e) in the hybrid expansion. Both calculations are 
technically feasible and should be done. 

In summary, we have analyzed the structure of the heavy-quark expansion for 
the inclusive, semileptonic B — > XI v decay rate with a lepton invariant mass cut 
<? 2 > <?o- This expansion is characterized by a hadronic scale /i c < m c determined by 
the value of g 2 . Because m b ^> /i c ^> A, the heavy-quark expansion can be organized 
as a combined (hybrid) expansion in two small mass ratios. The physics associ- 
ated with the two large scales mj, and /i c is disentangled using the HQET, whereas 
the physics on the scale \i c can be separated from long-distance physics associated 
with A utilizing an operator product expansion. We have used this formalism to 
obtain a RG-improved expression for the leading short- distance coefficient in the 
heavy-quark expansion at NLO. The summation of large logarithms in the hybrid 

2 We stress, however, that using the so-called Upsilon mass defined as one half of the mass of the 
T(IS') bottomonium state [31] does not eliminate the uncertainty associated with the variation of the 
6-quark mass. As discussed in [33], this choice obscures the presence of an unknown nonperturbative 
contribution to the bound-state mass, which is neglected in the perturbative expression of the B- 
meson decay rate in terms of the Upsilon mass. In other words, in such a scheme the value of m^, 
is known (by definition) with very high precision, but for consistency the uncertainty shown in the 
third column in Table 1 must then be added to the other theoretical uncertainties. 
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expansion turns out to be important and strongly enhances the overall size of the 
perturbative correction. We have also emphasized that in order to obtain a stable 
perturbative prediction it is important to eliminate the 6-quark pole mass in favor of 
a low-scale subtracted quark mass, such as the PS mass. Finally, we have presented 
several independent estimates of higher-order power corrections in the heavy-quark 
expansion, which at present do not permit a rigorous treatment. We find that with 
realistic values of the lepton invariant mass cut the overall theoretical uncertainty in 
the extraction of \V u b\ is about 10%, which is larger than previously estimated but 
still significantly less than the current uncertainty in this parameter. 
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